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Abstract. A Regge-inspired model is used to discuss the hard exclusive two-body hadronic reactions 
{pp tv'tv+ , tv^tt" , K'K+ , 'K°K") for the FAIR facility project at GSI with the Panda detector. The 
comparison between the differential cross sections predictions and the available data is shown to determine 
the values of the few parameters of the model. 
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1 Introduction 

The complex angular momentum method introduced by 
Regge [1] has been copiously used in particle physics. The 
Regge approach makes a non-trivial connection between 
particles, resonance spectra and high energy scattering 
amplitude behaviour. In the early 1970's, there was con- 
siderable interest in two-body hadronic reactions, both 
from a theoretical [2,3,4,5] and an experimental [6,7,8] 
view point. An important compilation of the data in the 
high energy regime were interpreted in terms of the Regge 
trajectories albeit with a limited succes. 
Later on, the Regge poles in connection with QCD [9] 
and the phenomenology based on the Pomeron hypothe- 
sis [10,11] led to a predictive approach in hadron-hadron 
scattering in the asymptotic high energy regime with a 
small number of parameters. 

The planned FAIR facility project at GSI together 
with the Panda detector [12] will allow one to reexam- 
ine the two-body hadronic reactions in antiproton-proton 
annihilation with unprecedented high accuracy. 

In principle, the different obscrvablcs for exclusive two- 
body hadronic reactions could conceivably be calculated 
from QCD, the fundamental quantum field theory of strong 
interaction. In practice, such obscrvablcs cannot be pre- 
dicted in a simple way with only quark and gluon degrees 
of freedom. A huge number of diagrams has to be calcu- 
lated as well as the contribution from pinch singularities 
which are related to the independent multiple scattering 
mechanism [3]. 

As said above, the Regge phenomenology is an alter- 
native and economical way to fit the available data of 



two-body hadronic reactions. The parametrization of the 
amplitudes in term of Regge trajectories is succesfuUy in- 
vestigated in photoproduction and clectroproduction of 
mesons and kaons [13,14,15]. 

It is well established that the Regge description is valid 
only for high energy (large s) but low negative momentum 
transfert {—t <^ s), where s, t and u are the usual Mandel- 
stam variables of the process. The linear trajectories are 
found to approach negative constants for large negative t 
[16,17] and the forward and backward Regge regimes must 
be joined smoothly onto the fixed large angle perturbative 
QCD behaviour. This extrapolation of the Regge model to 
the larger momentum transfer region is based on the sat- 
urating Regge trajectories {a{t) — > — 1 when i — > — oo ) 
and the introduction of the hadronic form factors in order 
to respect the dimensional counting rules [18,19]. 

In this paper, we investigate two pseudoscalar meson 
production in antiproton-proton annihilation. A phenome- 
nological model based on the Regge trajectories exchanged 
in the i-channcl and the vector meson in the s-channel 
is developed. An elegant procedure to connect the Regge 
trajectories in the two ranges of low and large negative 
momentum transfer t is discussed in detail and the com- 
paraison with the available data for two pseudoscalar me- 
son production is done. 



2 Model for two body hadronic reactions 

An attentive inspection of the two-body hadronic reac- 
tions in antiproton-proton annihilation data at > 3 
GeV shows a lack of theoretical investigations in both 
charged and neutral pseudoscalar meson pair production 
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[7,8,20]. Most of the existing studies for < 2 GcV use 
the effective lagrangien formahsm with the nucleon and 
Z\-resonance exchange picture. 

In the present paper, we perform a global analysis of the 
world data of two pseudoscalar meson production, using 
the Regge formalism in the high energy regime at invari- 
ant energy above 3 GeV. Wc reggcize the nucleon and A 
exchanges by appropriate Regge trajectories a{t) (Figs. 
1-2). 



2.1 Regge amplitudes 

To illustrate our formalism in its technical detail, we write 
down explicitly the amplitudes of the two meson produc- 
tion in antiproton proton annihilation. The vertex struc- 
ture of theses amplitudes is kept in the usual form given 
by the Feynman diagrams, using the phenomenological 
Lagrangians ( in appendix , restricted to energies up 
to the octet-baryons excitation region. Such a procedure 
was adopted in [13,14] for photo and electroproduction of 
mesons. 



TT+ 



The amplitude of the nucleon-exchange diagram in the 
f-channel reads : 



, 77-+ 




Fig. 1. Feynman diagrams for pp ^ -k tt"* 



where the vertex function C'^^^ contains the coupling 
constant and the form factor. The fermion propagator is 
the usual Feynman propagator : 



„2 _ ]\p 

Qn is the four-momentum of the exchanged nucleon. 



And for the zi-exchangc diagrams in the t-channcl: 



A good description of the forward and backward be- 
haviour of the differential cross sections of two-body hadronic 
channels was achieved in a satisfying way. Our original 
idea is to perform also a systematic analysis of the data 
in the large angle scattering region at large negative mo- 
mentum transfer t, which is the domain of predilection of 
perturbative QCD. To accomplish this delicate work, we 
assume the saturating trajectories for t — > — oo and the 
vector meson exchanged in the s-channel without reggeiza- 
tion to fit the available data. We only assume that Regge 
poles will control both the small |t| region and the large 
|t| region. Futhermore the Reggeon amplitude will domi- 
nate the Pomeron amplitude at fixed large angle and not 
too high energy. The details of this approach is described 
below. 



AA{s,t) = vp C^^^ V^^iQA) C:^^ Up (2) 
The propagator of the A is the Rarita Schwinger one : 

/lA + Ma 



V'^iqA) 
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The isospin coefficients in the amplitude of each diagrams 
in Figs. 1, 2 arc displayed in table 1. 



Table 1. Isospin coefficients in the amplitude of each diagram. 
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Fig. 2. Feynman diagrams for pp — >■ tt^tt" 



The basic idea of the Regge approach formalism is to 
replace the particle exchange Feynman diagram with a 
mass m by the exchange of certain quantum numbers. 
In the channels of our interest, this corresponds to the 
exchange of dominant baryon trajectories in the t and u 
channel and the usual Feynman propagator l/(t — m?) 
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where m. the mass of the exchange particle, is replaced 
by the Regge propagator. 



1 



t — 



1 + 5 exp{—i'!Ta(t)) 
sm{TTa{t))r{a{t) + 1) \so 



(3) 



where sq ~ 1 GeV^, 5 = ±1 is the signature of the tra- 
jectory and a{t) — aQ + a't is the Regge trajectory. The 
Regge propagator of eq. (1) is reduced to the Feynman 
propagator if t — )■ m^. It must be noted that formula (3) 
is only valid for a non-degenerate trajectory of spinless 
particles. 

The main ingredient in the Regge model is the trajec- 
tories defined by the spin J and mass mj of the particles 
with a fixed G-parity and the relation a{mj) = J. For 
two pion production in antiproton proton annihilation, 
the trajectories pertinent to our approach are shown in 
Fig. 3 . It is well known that Regge trajectories can be 
either non-degenerate or degenerate. According to Figs. 1 
and 2, the scheme adopted in our approach for two pion 
production in antiproton proton annihilation is discussed 
below. 

In order to fit the available data, we assume a degen- 
erate trajectory for the nucleon exchange with ajv(^) = 
—0.37 -I- 0.98 t and the Regge propagator 



exp(-i7r(ajv(i) + 1/2)) 



sin(7r(a7v(i) + 1/2)) r(aw(t) + 1/2) \so 



aiv(t)-l/2 



(4) 



For the Z\-exchange trajectory, a non-degenerate tra- 
jectory is used in both t and u channels a/i(i) = 0.1+0.93 1 
as follows 



1 — cxp{~in{aA{t) — 1/2)) 7ra^ /s 

2sm{TT{aAit) - 1/2)) r{aA{t) - 1/2) \^ 



2.2 The model 



azl(t)-3/2 

(5) 



For the large —t and —u regions, the linear Regge trajecto- 
ries fail to reproduce the data. This is the hard scattering 
region where the reaction mechanism exhibit a different 
behaviour. The differential cross sections are more or less 
independent of t and exhibit a scaling behavior with en- 
ergy as predicted by the dimensional counting rules. When 
extrapolated towards the large angle region scattering, one 
needs to introduce a new parameter tsat, the nominal value 
of t where the saturing trajectory {a{t) — >■ —1) becomes ef- 
fective. The original procedure of this paper is the smooth 
connection between the soft and hard mechanisms in or- 
der to always respect the concept of Regge trajectories 
(see fig. 3). 

Now the question can be raised how to calculate the 
parameter a' of the linear Regge trajectory in the transi- 
tion region between the hard and soft regions. 
Normally, a' is the slope of the linear Regge trajectory 



-0.5 - 



-1 

t (GeV/c)' 



Fig. 3. Nucleon Regge trajectory of the model in the full kine- 
niatical range of t 



ill the t > region. In the intermediate region between 
the asymptotic large t < tgat and t < 0, we assume the t- 
dependency of a' is determined by the slope of the straight 
line connecting the corresponding point on the trajectory 
and those with (a(m^) = J). 

The Regge amplitude is factorized in terms of a Regge 
propagator and a vertex function similar to the particle 
exchange Feynman diagram. The vertex function contains 
the form factors of the exchanged particles. The monopole 
form is adopted in the hard scattering regime where the 
form factors of the outgoing particles are crucial to re- 
cover the counting rules and the dipole forms for baryon 
exchanged in the soft region. A free parameter |xo| is in- 
troduced to separate the two kinematical regions. 

For I a; I > |xo|, we introduce a monopole pion form factor 
which corresponds to the hard scattering regime 



NBti 



B = Nov A 



with yl^^^ = 0.85 GcV, Aj^^^ = 0.55 GcV and 

For \x\ < \xq\ a dipole form factor is introduced both for 
the nucleon or the A 



A' -M 

B a 

Al-x 



2 \ 2 



The value of the parameter A^ is fixed by the constraint 
dX x = xq 

F^{xo)=Fl{xo) 

For two charged pions channel, the p-exchange dia- 
gram in the s channel is investigated with the traditional 
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Feynman formahsm to write down the corresponding am- 
phtude. Guided by the previous works [13,14], we do not 
reggeize the p exchange. The vertex functions of p-K-K and 
pNN are given by the standard effective Lagrangian with 
the coupling constants and form factors assumptions : 



L = —K 



7, 



(6) 



The scalar product in isospin space is explicitly written 
as: 

T p" ^ ^/2 [>^p /3+'' + Wn p *'p] 

+'^p P°'' ~'^n p'''' (7) 

The pTTTT Lagrangian is taken from ref. [22] 

^P^^ = 9p^^ [5" 71- X7v]-p^ (8) 
which gives in terms of the physical pion fields 

+ (9"^° TT--^"*- n°)p+ } (9) 



with the values of the coupling constants: 

^ V4^* 0.317, Kp = 6.033, = ^4^^ 
The usual p propagator without reggeization is: 



P'^p'/M^p 



s - M'p + iMpFp 



determine the values of some of the free parameters of the 
model. Let us emphasize the complimentary data in the 
charged and neutral pion channels. At high energy, one 
has only the data in the large angle scattering region for 
neutral pions instead of forward and backward regions in 
the charged channel. This leads to a strong constraints for 
fitting the data. 

As can be seen in the figures 4 and 5, the main features 
of the angular distribution are reasonably well reproduced 
as well as for the neutral pions as for the charged pions. 



The product of the form factors in the s-channel reads 



with ApT^TT = 0.95 GeV. is a normalization constant 
close to 1 taken as a free parameter. The value adopted in 
our calculation is = 0.9. 

To close this section on the model, we would like to em- 
phasize the small number of the free parameters (only 7) 
with our Regge approach to fit the two pion channels. 
These parameters are : tsat.N, isat,A for saturing trajec- 
tory; ^jv7Vir,^JV7V/i for the form factors; xqtv = for the 
transition of the dipole to the monopole form factor of the 
exchanged particule and the outgoing particle and finally 
ApTTTT, A^ for the p diagram contribution. We will deter- 
mine the values of these parameters in the next section 
after a best fit with the available data. 

3 Cross section predictions and comparison 
with available data 

Our predictions of the two pseudoscalar meson produc- 
tions can be compared to the available data, in order to 
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-1 -0.8 -0,6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 

cos 0c„ 

Fig. 4. Differential cross sections for pp — >■ tt^tv^ . Data from 
[7,8,21], red dashed line for s — 11.3 GeV^ and black solid line 
for s = 13.54 GeV^ 
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In this latter case, the asymmetry in the experimen- 
tal cross sections between forward (F)) (cos(0) ~ 0.90) 
and backward (B) (cos(6') ^ —0.90) angles are different 
at s = 11.3 GcV^ ( F/B- 3/2) and at s = 13.54 GcV^ 
( F/B^ 1/0-2) and this difference as a function of the 
incident energy is not reproduced by the calculation. At 
s = 11.3G'eT^^, a shoulder in the angular distribution at 
cos{9) ^ 0.65, which is missing at higher energy, is also 
not reproduced by our calculation. 

In Figs. 6 and 7, we display the contributions of each 
diagram of the two pion channels. 
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Fig. 6. Contribution of different diagrams in Fig. 1 for pp — > 
■k~tt'^ at s = 13.54 GeV'^ . Red line : nucleon, blue line ■.A'^'^ , 
green line :A~^ and pink line : p contributions 



In the case of charged pions, the forward scattering is 
dominated by the nucleon trajectory exchanged and the 
main contribution of the backward scattering is the Z\"*"+ 
in the it-channel for tt^tt"*". This is the reason of a dissym- 
metry only observed in the cross section for two charged 
pions production. Let us emphasize that the cross section 
ratio between the Z\++ contribution at cos 9 = —1 and the 

at cos 9 = +1 in Fig. 6 is given by the square of the 
isospin coefficients in Table 1. The contribution of the p in 
the s-channel is small, but its interference with the others 
contributions could play a role to fit the data in the large 
angle scattering region. 



In the case of neutral pions, the scattering is domi- 
nated by the nucleon trajectory for |cos0| < 0.65, while 
for I cos9\ > 0.80 the A and nucleon contributions are the 
same order of magnitude. In contrast with the charged 
channel, this feature could be checked with forthcoming 
data at forward angle. In Fig. 8, we display the compar- 
ison between the cross section prediction of our model 
versus the quark interchange model with a normalisation 
taken from the cross channel namely the pion-proton elas- 
tic scattering [8] . The main difference is the finite value of 
the cross section at cos 9 ~ ±1 and the large angle scatter- 
ing —0.6 < cos 9 < 0.6 where our model predicts a larger 
value of the cross section. Indeed, the quark interchange 
model is not appropriate to fit the data at forward and 
backward scattering and the lack of data near the 6* = 90 
degrees does not allow one to discriminate between the 
two models. 
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cos 0^.„ 

Fig. 8. Regge approach versus Quark interchange model [2] 
for pp — >■ 7r~7r"'" at s = 13.54 GeV^ . Red dashed line : Quark 
interchange model; black solid line : Regge approach model. 



In summary of this section, we display in the table 
2, the values of different parameters of the model. These 
values are obtained with a reasonable fit without stan- 
dard minimization procedure, to leave room for future 
constraints from data. 



Table 2. Values of different parameters defined in the text for 
two pion channel 



tsat.N 


tsat.A 




An An 


Xo,N — Xo.A 






-2.6 


-3 


0.85 


0.55 


-1.4 


0.9 


0.95 


GeV^ 


GeV^ 


GeV 


GeV 


GeV^ 




GcV 



4 Two kaon production 



With the Regge approach formalism, the two kaon produc- 
tion channels are dominated by the exchange mechanism 
of the baryons A and S in the i-channel and the ^-meson 
in the s-channel. We display in Fig. 9 the corresponding 
Feynman diagrams. 

In the case of strange meson production, the exchanged 
baryons have to carry quantum numbers at each vertex in 
order to respect the strangeness conservation. 




01 23456789 10 

-t {CeV/cf 



Fig. 10. Differential cross sections for pp — >■ K K . Data 
from [7,8], red dashed line for s = 11.3 GeV^ and black solid 
line : s = 13.54 GeV^ 



1.° 




P 



Fig. 9. Feynman diagrams for pp ^ K K 



The Lagrangians used in our model, the coupling con- 
stants and the form factors are given in the appendix VK\ 
Assuming the same set of free parameters of the Model 
given in table 1, plus three new parameters namely 
A-NAK = Anek — 1-05 GeV for the form factors; and 
finally A^kk = 1.3 GeV, = 1 for the <!> diagram con- 
tribution , we are able to fit reasonably the charged kaon 
channels. As in the case of charged pions production, at 
s = 11.3 GeV^, the main discrepancy between the data 
and our caculation is located at t ~ —2 GeV^ correspond- 
ing to cos(6') ^ 0.65. The agreement is better at higher 
energy (s = 13.54 GeV^). 

In contrast to the two pion channels, we have no contribu- 
tion from the u-exchanged diagram. This is an important 
check of the model. It seems that the data confirms this 
fact, but the error bars are too large to make a conclusion. 
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-t (GeV/c)' 

Fig. 11. Contribution of different diagrams in Fig.9 for pp — >■ 
K~ for s — 13.54 GeV^. Green line : A contribution; blue 
line : S'^ contribution and pink line: ^ contribution. 




Fig. 12. Feynman diagrams for pp — >■ K K 



5 Conclusions 

In our model, there is a nice continuity between the soft 
and hard momentum transfer regimes. Both are controlled 
by the Regge pole trajectory exchanges and the saturing 
trajectory at large —t. With few free parameters ( 10 for 
two pseudoscalar meson production) and the Regge ap- 
proach, we are able to fit reasonably well the data in an- 
tiproton proton annihilation. 

It seems that our calculation reproduces correctly the data 
for s > 13 GeV^ and with less succes at a lower inci- 
dent energy (s — 11.3 GeV^). Nevertheless, we have to 
make the following remarks: i) the discrepancy with the 
data measured with the same experimental set-up occurs 
mainly at the same angular range for both charged pions 
and kaons. ii) the experimental error bars in the publi- 
cations show large uncertainties in the data, iii) In the 
case of the charged pions, the variation with the incident 
energy of the Forward-Backward cross sections should be 
confirmed experimentally. Therefore, we need more data 
in the whole angular range and in a large kinetic energy 
domain with better statistics and small systematic errors. 
Large angle scattering data will provide constraints on the 
parameters of the model. In the hard scattering region, 
reliable perturbative QCD calculation of these channels 
would be very valuable towards a consistent conclusion. 

The Panda detector gives an unique opportunity to 
measure these cross sections with very small statistic er- 
rors in an energy range 5.4 GeV^ < s < 27 GeV^. For 
pp — >■ 7r~7r+, the expected PANDA counting rates for one 
week of beam time at full luminosity L — 2.10~^^ cm~^ 
s~^ within a 0.1 wide cos6' bin at = 60 degrees are dis- 
played in table 3. 



T 




10 r 



2 4 6 8 10 

-t (GeV/c)' 

Fig. 13. Differential cross sections for pp — > K^K'^ with each 
diagram contribution of Fig. 12 at s = 13.54 GeV^. Blue line : 
contribution; pink line : contribution. 



Table 3. Expected counting rates of pp — tt tt"*" for different 
values of s. 



s (GeV^) 


13.5 


16 


20 


Event numbers 


10^ 


2 10^ 


4 10* 
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A Lagrangians 

We give the Lagrangians used in the t or u-channels with 
the appropriated coupling constants. The adopted conven- 
tion is that the charge associated with the field corre- 



8 



J. Van de Wiele et al.: Regge description of two pseudoscalar meson production in antiproton-proton annihilation 



sponds to the annihilation of the corresponding incoming 
particle. 



A.l Pion production 

We have used the pseudo- vector ttNN lagrangian which 
is the one most commonly found in the literature : 

C^NM = -KIn^nI'^^ ^ ■ l^"^) (A-1) 



Or cxplicitely, in terms of charged pion fields : 

pv pv 
r — _J^^ y 

irNN ttNN 

{72 [^p 7^7^ (9" 7i+) W„ + W„ 7%^ (9" 

+ t'Im (9" tt") - 7^7^ (a" 7r°)<I'n } 

(A-2) 

with the coupling constants [26] 



Knn = ff^««/2A'/„, = * 12.562. 

The ttNA Lagrangian is given in a compact form sim- 
ilar to cq. I A- II : 



The coupling constants gP*^ . and q'^^ arc much less 
constrained experimentaly than g^j^ ■ The algebra asso- 
ciated with the two octet representations of SU(3) gives a 
link between these coupling constants [23]. We have cho- 
sen the values given in [13,27] with the upper and lower 
limits corresponding to 20% breaking of the SU(3) sym- 
metry : 



ps 
NA 
ps 

NAK min 
ps 

NEK 
ps 



0^° = \/47r (qP^ 
^NAK ^ NAK min 



■9 



ps 

NAK max 



)/2, 



-9/2, g 



NAK max 
ps 



47r gP*^ + q 

NEK min ^ NEK max 

0.9, = 1.3. 



)/2 



NEK min ' NEK max 

The NNip Lagrangian is given in [24] 



with K 



= 0.5. 



(A- 



C 



NN4. = -0.90, rv0 
The KK<P Lagrangian is given in [25] 



(pKK 



V2 



K+ ff^K- - R-d^K+ 



(A-9) 



^.NA = Kna K ( ■ ) + h.c (A-3) 
with [26] K^^^ = g^^jM^, g^^^ = 2.13 



r 6.46 for K- K+ 
~ 1 6.67 forifo K° 



A. 2 Kaon production 

Following previous calculations, we assume a pseudo-scalar 
form for the baryon-baryon-meson lagrangian [23] 



NAK ~ ^NAK 



+ ¥n 7' A'" + l^^"^p } (A-4) 



The NEK Lagrangian is written as a sum of two terms : 



L 



= gP^ 

NEK ^NEK 



^NEKi') + ^NEK('') 



with 



V2[S ■y^'R-'Pn + E 7^7^^^] 
+ ^\'' k- % ¥° >P„ (A-6) 



and 



£l"^^ {II) = V2[>Pp-f"K°E+ + 7"^ K+ E- 



+ ifp 7^ k+ e" - 7^ k° r° (A-7) 
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